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MEASUREMENT  OF  2D  WAVEFRONT  TILT 
USING  ID  ARRAY  DETECTORS 

Jiang  Zhiping,  Zhong  Hairong, 

Lu  Qisheng,  and  Liu  Zejin 

Department  of  Applied  Physics 

National  University  of  Defense  Technology 

Changsha  410073 

Abstract:  A  method  that  uses  a  one-dimensional  array 

detector  to  measure  the  two-dimensional  wavefront  tilt 
of  a  beam,  is  proposed.  A  preliminary  experiment  is 
given  to  verify  the  principle. 

KEY  WORDS:  wavefront  tilt,  array  detector. 

1.  Introduction 

The  wavefront  tilt  of  a  beam  is  required  to  be  measured  in 
numerous  applications,  such  as  adaptive  optics,  beam  propagation 
through  the  atmosphere,  automatic  collimation  of  beams  and  laser 
beam-based  measurement  of  the  tilt  of  objects,  etc.  Generally 
speaking,  all  these  problems  are  two-dimensional  problems,  in  which 
the  measurement  of  two-dimensional  wavefront  tilt  is  required. 

Tilt  measurement  can  be  equivalent  to  the  measurement  of  the 
beam  intensity  centroid,  which,  in  earlier  times,  was  realized 
using  four-quadrant  detectors.  With  the  development  of  array 
devices  (such  as  CCD  photoelectric  detectors),  however,  the  four- 
quadrant  detectors  are  now  replaced  by  array  detectors. 
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To  measure  the  two-dimensional  light  intensity  centroid,  a 
face  array  device  is  required,  which  consists  of  numerous  detectors 

I 

(e.g.,  a  face  array  CCD  contains  hundreds  of  thousands  of 
photosensitive  elements).  In  such  devices,  data  are  transmitted 
sequentially  and  therefore,  the  frame  frequency  cannot  be  high, 
such  as  25Hz  in  the  CCD  camera,  which  is  insufficient  for  rapid 
measurement. 

Compared  with  face  array  devices,  the  linear  array  devices 
involve  many  fewer  detectors  and  have  higher  a  frame  frequency. 
Now  the  problem  is:  Is  it  possible  to  use  linear  array  devices  for 
measuring  the  two-dimensional  wavefront  tilt?  A  technique  of 
measuring  the  two-dimensional  wavefront  tilt  using  an  one¬ 
dimensional  array  detector  is  proposed  in  this  paper;  its 
feasibility  was  verified  through  a  preliminary  experiment. 

2.  Principle 


A.  W.  Lohmann  et  al.  pointed  out[l]  that  by  placing  a 
specially-distributed  curved  aperture  at  the  incident  plane  of  a 
2f  optics  system  (front  focal  plane) ,  two  one-dimensional 
functions,  namely  the  Airy  function  and  the  Daguerre  polynomial  can 
be  derived  at  the  output  plane  {  rear  focal  plane)  ,  With  the 
instructions  from  that  paper,  we  were  considering  a  problem,  that 
is:  Can  we  use  Lohmann 's  technique  in  measuring  the  wavefront  tilt 
of  a  beam?  Fortunately,  the  two-dimensional  wavefront  tilt  can  be 
changed  into  one— dimensional  as  long  as  a  dual  straight  line 
aperture  is  employed. 
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Fig.  1.  Schematic  diagram  of  the  optical  setup 

Take  the  2f  optics  system  shown  in  Fig.  1  as  example.  Let  the 
incident  light  be  a  plane  wave  that  can  be  expressed  as 
\  cxp[iy(0.i +  ,  where  0^  and  Oj ,  respectively,  are  helical 

angles  along  x  and  y  directions.  Suppose  the  transmittance 
function  of  the  input  plane  is  T(x,y},  then  the  optical  amplitude 
of  the  output  plane  E(5,ii)  can  be  calculated  by  using  the  Fourier 
transform  as  follows: 

Jr(i,jr)exp[iy(0,j:  4- 0,y)]exp(» 

where  X  is  wavelength;  f  is  the  focal  length  of  lens;  a  is  the 
dimension  of  lens.  Here,  for  convenience,  let  the  lens  be  square, 
which  in  fact  is  not  square,  and  let  transmittance  be  dual  straight 
lines,  then 

T(i, j)  =  d(y  —  Cl*)  +  ~  ‘^i*)  { 2 ) 

where  the  two  constants  c^  and  Cj,  respectively,  are  the  slope  of 
the  two  straight  lines.  By  introducing  Eq.  (2)  in  Eq.  (1)  ,  the 
following  can  be  derived: 

Eii,,,) _ Jj  J  j  Cd(?  -  c,*)  +  Si9  -  c**)>xp(i  f  +  0,9  +  ^^^))dxd, 

— «  — « 

=  j  cxp(iC.x)dx  j  exp(iC^)d*)  =  ^Csinc((7,<i)  +  sinc((7:<i)D  ^  ^  ^ 

where 

+  +  j_,,2  (4) 

From  the  nature  of  the  sine  function  it  is  known  that  in  Eq.  (3) , 
if  a-^,  then  sinc(«i) -►  i(*)  ,  so 

F(4,7)1,—  =  ^CdCC.)  +  HCtB 
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(5) 


In  Eq.  (5),  the  definition  of  function  3(2)  is  3(0)  —  lf3<*  ^  0)  0 

The  non-zero  E{C»n)  iri  Eq.  (5)  should  satisfy  the  following 

i 

expression: 

<7i  =  0,  » =*  1,2  (g) 

At  the  output  plane,  the  intensity  distribution  is  represented  by 
two  bright  straight  lines.  From  Eqs.  (6)  and  (4) ,  the  two  straight 
lines  can  satisfy  the  following  formula: 

+  +  .=  1,2  (7) 

and  (5^,iij^)are,  respectively,  the  points  on  the  two  straight 
lines  at  the  output  plane.  By  solving  Eq.  (7) ,  the  following  can 
be  derived: 

_  1  »?|Ct  —  +  C|Ci(i?i  —  ih) 

_1_  —  ^1  +  —  7iCi 

'  /  C, -c. 


Eq.  (8)  implies  that  the  two-dimensional  wavefront  tilt  can  be 
calculated  if  only  the  coordinates  of  the  points  on  the  two 
straight  lines  are  measured,  which  can  be  realized  with  an  one¬ 
dimensional  array  detector.  For  the  convenience  in  use,  Eq.  (8) 
can  be  simplified.  First,  let  a  linear  array  detector  be  parallel 
to  the  ^-axis,  then  Hj=ii2=constant .  On  the  other  hand,  it  is  not 
difficult  to  understand  that  the  measurement  of  wavefront  tilt 
focuses  on  its  variation;  therefore,  taking  the  dynamic  tilt  into 
account,  Eq.  (8)  can  be  written  as: 


X(0  =  4 


f  Cl—  Cj 

J_  .^,(0  -  .^.(O 

/  Cj  —  Cl 


where  t  is  time  variable;  all  the  quantities  associated  with  time 
in  Eq.  (9)  refer  to  the  variation  relative  to  the  time  t=0.  For 
instance,  -^.(O  =  d.(0  —  ff.(O)  .  In  addition,  if  C2= — Cjsl  is 

evaluated,  then 
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(10) 


-^^>.(0  =-  ^CA4,(0  + 

-^,(0  =  -  ^*C0D 

The  linear  array  detector  can  also  be  placed  parallel  to  axis  n  . 
In  this  case,  the  formula  derived  resembles  Eq.  (10)  but 
essentially  they  are  identical. 

In  actual  measurements,  a  cannot  be  infinite,  and  Eq.  (3) 
cannot  be  simplified  as  Eq.  (5),  yet  the  function  sinc(x)  has  a 
main  maximum  at  x=0.  It  should  be  understood,  therefore,  that  the 
foregoing  discussion  is  directed  at  the  major  maximum. 
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Fig.  2.  Intensity 
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distribution  at  the  output  plane 

Fig.  2  is  a  schematic  diagram  showing  the  intensity 
distribution  at  the  output  plane  and  the  location  of  a  linear  array 
CCD  (imaginary  line  in  the  figure)  .  Since  the  wavefront  of  the 
incident  light  is  generally  tilted,  the  two  bright  lines  at  the 
output  plane  do  not  cross  at  the  original  point. 

3.  Experimental  Verification 

To  verify  the  feasibility  of  the  foregoing  concept,  we 
conducted  a  simple  and  preliminary  experiment  to  measure  the  wedge 
angle  of  an  optical  wedge.  The  experimental  setup  is  shown  in 
Fig.  1.  The  aperture  in  Fig.  2  was  realized  through  photography; 
through  photography,  the  two  0.3-mm-wide  mutually  vertical  straight 
lines  were  minimized  by  a  factor  of  10,  and  the  aperture  thus 
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derived  had  a  width  O.OSmm  (theoretically,  the  smaller  the  width  of 
aperture,  the  better;  however,  considering  its  effect  on  the  light 
intensity  and  for  the  convenience  in  making  the  aperture,  it  is 
appropriate  to  select  it  as  approximately  0.03mm). 

The  laser  beam  emitted  with  an  He-Ne  laser  device,  propagating 
approximately  10  m,  arrived  at  the  narrow-slit  aperture  located  at 
the  input  plane  (front  focal  plane  of  the  lens) .  The  beam  radius 
at  the  aperture  was  approximately  10mm,  and  the  focal  length  of  the 
lens  was  approximately  1.08m.  A  visible  light  linear  array  CCD  was 
employed,  with  an  interval  of  photosensitive  elements  14Mm,  and  a 
number  of  image  elements  2048.  The  data,  through  an  A/D  converter, 
were  processed  on  the  computer. 


Position 


Fig.  3.  CCD-detected  intensity  distributions 
(a)  original  -  (b)  a  glass  plate  is  inserted 
(c)  the  glass  is  rotated  by  approximately  90^* 

Fig.  3  shows  the  intensity  distribution  detected  with  a  CCD, 
the  abscissa  being  the  location  of  CCD  image  elements.  Fig.  3(a) 
shows  the  intensity  distribution  without  glass  plates,  while 
Fig.  3(b)  shows  the  intensity  distribution  with  a  glass  plate 
inserted  in  front  of  the  aperture  at  the  input  plane;  then,  by 
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rotating  the  glass  plate  by  90“,  the  intensity  distribution  was 
measured  again  [Fig.  3{c)l.  In  calculating  the  centroid  of 
intensity  distribution,  the  intensity  distribution  was  properly 
treated  with  thresholds.  The  glass  plate  has  a  dip,  just  like  the 
tilt  of  incident  light. 

Compared  with  Fig.  3(a),  the  intensity  centroids  of  the  two 
bright  lines  in  Fig.  3(b)  are  close  to  each  other,  with  much  the 
same  travel  motion,  and  it  is  known  from  Eq.  (10)  that  the 
wavefront  tilt  primarily  goes  along  the  y  direction  at  this 
instant.  The  intensity  centroids  of  the  two  bright  lines  in 
Fig.  3(c)  simultaneously  move  to  the  left  indicating  that  at  this 
instant,  the  tilt  mainly  goes  along  the  x  direction  [Notice:  the 
glass  plate  in  Fig.  3(c)  has  been  rotated  by  approximately  90® 
relative  to  Fig.  3(b),  which  conforms  to  the  measurements  here]. 
The  wavefront  dips  derived  after  data  processing,  respectively. 


are : 

^0.  =0.00' 

=1.92' 

corresponding 

to 

Fig.  3(b) 

(11a) 

=1.99' 

=0.00' 

corresponding 

to 

Fig.  3(c) 

(11b) 

The  dip  of  the 

glass  plate 

is 

approximately  4.0', 

and  its 

reflectivity  is  approximately  1.5;  thus,  the  introduced  dip  is 
(1 . 5-1) x4 . 0 ’ =2 . 0 ' ,  basically  in  agreement  with  Eq.  (11).  This 
result  shows  that  our  measurement  is  reliable,  and  our  technique  is 
feasible.  When  0«1,  we  have  the  expression 

During  the  experiment,  when  the  wavefront  tilt  introduced  by 
the  optical  wedge  was  not  exactly  along  either  the  x  or  y 
direction,  the  measurements  conformed  to  Eq.  (12)  .  Since  this  was 
a  principle-based  experiment,  we  did  not  take  vigorous  approaches 
to  eliminate  various  errors.  The  major  error  sources  include:  (1) 
The  glass  plate  was  not  exactly  perpendicular  to  the  optical  path. 
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and  the  wavefront  tilts  introduced  by  the  glass  plate  in  two 
experiments  Fig.  3(b)  and  (c)  were  different,  which  served  as  the 
major  reason  for  the  errors  in  Eqs.  (11a)  and  (11b);  (2)  the 
aperture  was  not  placed  at  exactly  45^^  (i.e.,  C2=-Cj=l  could  not  be 
strictly  ensured) .  This  is  a  systematic  error,  which  can  be 
eliminated  in  actual  applications;  (3)  since  the  entire 
experimental  setup  was  placed  on  an  ordinary  disk,  the  mechanical 
vibration  system  itself  carries  a  dynamic  wavefront  tilt. 

It  is  to  be  noted  that  this  technique  was  designed  to  measure 
the  wavefront  tilt,  but  it  measured  even  the  wavefront  tilt  caused 
by  mechanical  vibration.  With  this  technique  to  measure  the  two- 
dimensional  wavefront  tilt,  we  realized  the  two-dimensional  stable 
control  of  beams  [2]  .  Further  studies  and  applications  are  now 
being  underway. 

4.  Conclusions 

This  paper  proposed  a  technique  for  measuring  the  two- 
dimensional  wavefront  tilt  using  linear  array  detectors,  and 
verified  its  feasibility  through  an  experiment.  This  experiment 
demonstrated  that  this  technique  features  high-speed  measurements. 
The  frame  frequency  of  linear  array  detectors  are  over  2  orders  of 
magnitude  higher  that  of  face  array  detectors;  accordingly,  their 
measurement  speed  also  can  be  more  than  two  orders  of  magnitude 
higher  (for  instance,  the  speed  can  reach  a  few  kHz) ,  which  makes 
these  detectors  available  for  cases,  in  which  fast  measurement  is 
required.  In  principle,  this  technique  can  provide  sensitivity  and 
accuracy  as  high  as  the  sensitivity  and  accuracy  offered  by  linear 
array  detectors  in  one-dimensional  measurements  as  indicated  in  Eq. 
(10) .  The  price  of  such  detectors  probably  is  not  important,  but 
they  can  be  used  as  substitutes  when  only  linear  array  devices  are 
available . 

However,  this  technique  also  suffers  from  some  disadvantages. 
For  instance,  it  may  cause  light  intensity  loss,  for  which  it  is 
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only  applicable  for  high-power  beams  (if  two  mutually  vertical 
narrow  silts  are  replaced  with  two  sets  of  mutually  vertical  narrow 

i 

silts  with  numerous  parallel  narrow  silts  in  either  set,  then  the 
utilization  rate  of  incident  light  can  be  increased.  Such  kind  of 
aperture  can  be  derived  through  photography) .  Although  this 
technique  cannot  completely  replace  two-dimensional  detectors, 
still  it  is  useful. 

The  analysis  and  experiment  presented  in  this  paper  are  all 
preliminary.  Many  problems,  including  the  effect  of  nonplanar 
wave,  are  expected  to  be  solved  in  the  future.  The  ideas  of 
Lohmann  et  al.  as  well  as  the  concepts  of  the  authors  are 
noteworthy.  What  else  can  we  do  with  other  apertures?  The  authors 
hope  that  this  paper  can  somewhat  serve  as  a  probing  attempt,  which 
will  invite  more  valuable  papers  to  come. 
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TEN-THOUSAND  WATT  LASER  LARGE  AREA  AUTOMATIC  SPREADING  SYSTEM 

Among  the  "Eighth  Five-Year  Plan"  key  projects,  the  laser 
machining  center  affiliated  with  Tianjing  Textile  College 
successfully  developed  an  automatic  10,000-W  watt  large-area  laser 
spreading  system,  which  consists  of  a  10  kw  laser-based  spreading 
wide-band  scan  rotating  mirror,  and  a  laser  spreading  automatic 
powder  spray  machine.  The  scanning  mirror  has  two  functions:  wide 
band  and  focusing  scan,  with  a  scanning  width  10-50  nm  continuously 
adjustable,  and  a  scanning  frequency  of  60 , OOOtimes/second 
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continuously  adjustable.  In  the  scanning  band,  the  beam 
distribution  is  uniform.  The  automatic  powder  spray  machine  adopts 
adjustable  rotating  drum  material  level  measurement,  combined  with 
non-carrier  gas  powder  delivery,  the  powder  delivery  being 
continuous,  uniform  and  stable.  The  powder  quantity  delivered  is 
large  and  can  reach  0.5-200g/min  continuously  adjustable;  the 
precision  of  powder  delivery  is  52%;  powder  utilization  rate  can 
reach  98%.  The  materials  to  be  delivered  cover  a  wide  range, 
including  not  only  ordinary  Ni  base,  Fe  base  and  Co  base,  but  also 
super-fine  ceramic  powder  and  super-heavy  WC  powder.  The  single¬ 
channel  spreading  width  {10~35nm)  and  single-channel  spreading 
thickness  (0.2-8  nm)  allow  a  wide  range  of  adjustment. 

With  the  foregoing  rotating  mirror  and  powder  spray  machine, 
a  long-time  laser  spreading  test  was  conducted  on  a  #500  large 
roller.  The  test  shows  that  the  above-mentioned  system  worked  over 
10  consecutive  hours  with  8-kw  CO2  laser  spreading,  and  its 
reliability  index  met  the  requirements  from  the  State  Reliability 
Center. 

This  achievement  proves  to  be  a  significant  invention  in 
aspects  such  as  10,000-W  laser  wideband  rotating  mirror  linear 
spot-scan  technology,  powder  delivery  mechanism,  wideband  non¬ 
carrier  gas  powder  delivery,  and  other  measures  needed  for 
industrial  production.  Also,  this  system  well  conforms  to 
conditions  prevailing  in  China.  Its  overall  level  can  keep  up  with 
the  international  advanced  level.  This  accomplishment  was 
certified  by  the  related  department  of  this  country  in  January  1996 
and  began  to  be  applied  in  laser  spreading  and  quenching  with 
respect  to  multiple  rollers  below  #500,  automobile  rocker  arms, 
convex  axles  and  automobile  oil  seal,  bringing  about  great  economic 
and  social  benefits. 
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LIMIT  OF  ANISOPLANATISM  TO  LASER-GUIDE-STARS 

Fan  Chengyu  and  Song  Zhengfeng 

Anhui  Institute  of  Optics  and 
Fine  Mechanics 
Chinese  Academy  of  Sciences 
Hefei  230031 

Abstract;  The  main  limit  for  a  laser-guide-star  adaptive- 
optics-system  is  the  focal  anisoplanatism.  Firstly,  based 
on  the  filter  function  of  the  atmospheric  turbulence  sensing 
geometry,  the  variance  of  the  residual  wave-front  distortion 
resulting  from  focal  anisoplanatism  was  derived,  then  the 
effective  diameter  dj  of  the  imaging  system  was  derived, 
dj  was  calculated  according  to  the  atmospheric  turbulence 
model  in  Kunming  and  Xinglong.  The  calculation  shows  that 
the  correction  of  the  piston  and  tilt  is  important  for  the 
adaptive  optics  system. 

KEY  WORDS;  laser-guide-star,  adaptive  optics,  focal 
anisoplanatism. 

1.  Introduction 

Placing  an  adaptive  optics  system  (AOS)  in  the  traditional 
ground-based  astronomical  telescope  proves  to  be  an  important 
means  for  upgrading  imaging  quality.  In  this  case,  the 
acquisition  of  a  beacon  source  in  the  AOS  is  a  key  problem;  the 
beacon  is  supposed  to  lie  within  the  isoplanatic  angle  of  the 
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target  to  be  measured.  It  is  well  understood  that  in  the  range 
of  visible  light  and  infrared  light,  the  isoplanatic  angle  is 
extremely  small,  generally  only  a  few  angular  seconds.  In  this 
range,  there  is  no  naturally-originated  guide  star  (NOGS) 
available.  To  solve  the  problem  of  guide  stars  applicable  in 
astronomical  observations,  Feinleib[l]  in  1982,  and  Foy  and 
Labeyrie[2]  in  1985,  respectively,  proposed  the  concept  of  using 
Rayleigh  scattering  or  sodium  resonant  scattering  as  artificial 
guide  stars  (AGS)  in  AOS,  and  later  made  quite  a  few  successful 
demonstrations . 

Since  the  height  of  the  AGS  is  limited,  the  spherical  waves 
from  the  AGS  and  the  plane  waves  from  celestial  bodies,  while 
passing  through  atmospheric  turbulence  and  arriving  at  receiving 
apertures  along  different  paths,  can  form  an  anisoplanatic 
region,  or  so-called  focal  anisoplanatism. 

This  paper  primarily  discusses  the  limits  of  focal 
anisoplanatism  over  the  AOS.  First,  the  authors  derived  the 
formula  for  the  effective  diameter  dg  as  a  measure  of  the  focal 
anisoplanatism.  Then,  based  on  the  calculations  of  the 
atmospheric  turbulence  models  in  Kunming  and  Xinglong,  the 
authors  demonstrated  that  dg  is  closely  associated  with  the 
distribution  of  atmospheric  turbulence,  and  that  the  correction 
of  the  phase  distortion  of  low-order  terms  is  significant  to  the 
AOS. 

2-  Basic  Formulas 

According  to  the  analysis  made  by  Fried [3],  for  a  circular 
aperture  receiving  system,  a  spherical  wave  from  beacon  source 
and  a  plane  wave  from  a  celestial  body  can  generate  a  variance  of 
residual  wavefront  distortion  due  to  focal  anisoplanatism  as 
shown  in  Fig.  1,  which  can  be  expressed  as: 
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(1) 

where  D  is  the  diameter  of  the  aperture  of  a  receiving  telescope; 
dji  is  the  effective  diameter  of  the  AOS. 


object  at  infinity 


sodium  guide  star 
/Z  — 90km 


Turbclcnt  layer 

Rayleigh  guide  star 
//^10-30km 


Fig.  1.  Diagram  of  artificial  guide  stars 


0 

Here,  the  main  requirement  is  to  calculate  E‘;  taking 
advantage  of  the  systematic  integration  average  of  the  phase 
structure  function  at  the  plane  of  the  receiving  aperture,  it 
will  be  fairly  easy  to  derive  the  general  expression  for 
calculating  the  variance  of  the  residual  wavefront  distortion 
resulting  from  the  focal  anisoplanatism[4]  as  follows: 


F*  =  0,  207i-S  Jc*(2)[  J/(A:)j(K-,2)dK’jdz  ^ 

where  k|j=2n/A,  X  is  beacon  wavelength;  the  integration  of  z 
extends  from  the  receiving  aperture  to  the  celestial  body  of 
interest;  the  integration  of  K  is  a  two-dimensional  space  at  the 
plane  of  the  receiving  aperture;  f(K)  is  the  atmospheric 
turbulence  spectrum,  which,  in  our  case,  is  the  normalized 
Kolmogorov  turbulence  spectrum 


(3) 

g{K,z)  is  the  aperture  filter  function  representing  the  sensing 
geometry  of  turbulence,  which  takes  different  forms  with 
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reference  to  different  problems. 


First,  we  consider  the  variance  of  the  residual  wavefront 
distortion  under  the  height  of  the  beacon.  In  accordance  with 
the  analysis  made  by  Sasiela[4],  when  the  piston  and  tilt  are 
retained,  the  filter  function  takes  on  the  following  form: 


g(.K,z)  = 


-  z)/2H^\ 
KD(.H-z)/2H 


(4) 


By  introducing  Eqs.  (3)  and  (4)  in  Eq.  (2),  and  integrating  the 
azimuth  angle  of  space  K,  the  distortion  variance  in  the  presence 
of  the  piston  and  tilt  can  be  derived  as  follows: 

ES  =  2.  606ii  J  dzCKz)  J  dA'A'-'/’j  1  -  -  ^  ^  ^ 

0  0 

Finally,  through  simplified  calculations,  the  distortion  variance 
expressed  with  part  of  the  turbulence  moments  can  be  obtained  as 
follows 


^  =  0-  5007i-5scc^(Z)/^)*/V/i(^) 

Here,  the  definition  for  some  of  the  m-order  turbulence  moments 
is:  for  the  lower  part  of  the  beacon 

M 

=  Jdis-CKi)  (7) 

0 

and  for  the  upper  part  of  the  beacon 


(8) 


Normally,  researchers  are  interested  in  the  situation  after 
the  piston  and  tilt  are  corrected.  Based  on  Eq.  (2)  and  the 
filter  functions  of  the  piston  and  tilt,  the  distortion  variance 
caused  by  the  piston  and  tilt  can  be  calculated  separately. 
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Moreover,  as  long  as  the  variance  of  the  piston,  or  the  variance 
of  the  piston  and  tilt  is  subtracted  from  Eq.  (6) ,  the  residual 
variance  can  be  calculated  after  the  piston,  or  the  piston  and 
tilt  are  corrected. 


The  filter  function  of  the  piston  and  tilt  can  be  written  as 
follows [4] 


ftr  ^  J.(0/2)  JX.KDiH-z)nH')Y 

^,^^2  KDiH-z)/2H  I  (gj 

where  subscript  z  represents  piston  (P)  and  tilt  (T)  ;  it 
represents  the  piston  when  v=l  and  represents  the  tilt  when  v=2. 
Through  substituting  Eq.  (9)  in  Eq.  (2),  and  integrating  the 
azimuth  angle  of  space  K,  the  following  can  be  derived: 


=  1 


H 

3031-i  J drCl(5)4i-* 


JXKDiH  -  z)/2H'iY 
KDCH  -  z)l2H  I 


(10) 


Then,  by  introducing  the  variable  x=KD/2,  y=(H-z)/H,  and 
inserting  an  additional  term  with  zero  as  a  sum,  we  obtain 

E\z=  {  i=C\{z)l(z'>  (11) 

*"0 

/(.*)  =  Jd«-‘^(  [/;(*)  -  ^]+  -  5^]-  2[j.(*V.C*y)/y  -  5^*)) 

The  additional  term  is  introduced  in  Eq.  (12)  so  that  the 
integration  can  converge  when  the  piston  is  calculated;  a=l  is 
selected  in  calculating  the  piston,  and  a=0  in  calculating  the 
tilt. 


According  to  the  Mellin  transformation  technique  developed 
by  Sasiela[5],  the  variance  of  piston  distortion  can  be 
calculated  as 

E\,=  0.  0833iJsecj>Z)*/W  ^ ^ ^ 

and  the  variance  of  tilt  distortion  is 
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^  r  =*  0.  368kiscc^IJ^'*Mi  q  4  j 

By  subtracting  Eq.  (13)  from  Eq.  (6) ,  the  variance  of  residual  ! 
distortion  after  the  piston  is  corrected  can  be  derived  as 

^  M  =  ^sec^ZJ'-^CO.  5007mL(H}/H^'*  -  0.  0833^*^  ^  ^  ^  ^ 

By  subtracting  Eqs.  (13)  and  (14)  from  Eq.  (6),  the  variance  of 
residual  distortion  after  the  piston  and  tilt  are  corrected  can 
be  obtained  as 


5007Mli(H:>/H^'*  -  0.  iSlZui 


For  the  effects  of  turbulence  in  the  upper  part  of  the  beacon, 
the  variance  of  residual  phase  distortion  can  be  derived  by  using 
the  Zernike  model  proposed  by  Noll  [6].  After  the  piston  is 
corrected,  the  following  can  be  obtained: 

1.0299  (£>/r/)»/» 

Here,  the  atmospheric  coherent  length  r(jt,  calculated  above  the 
beacon  height,  can  be  expressed  as 


Hence , 


=  CO.  423ii^(^);)->/s 


(18) 


Similarly , 
corrected 


M  =  0.  4Z6kiseci^[fi'*ftj  (ff) 


(19) 

the  residual  variance  after  the  piston  and  tilt  are 
is 


i.e. , 


(20) 


^  =  0.  057iSsecj>Z>*'V  (H)  (21) 

Thus,  by  adding  Eqs.  (15)  and  (19),  the  distortion  variance 
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caused  by  focal  anisoplanatism  after  the  piston  and  tilt  are 
corrected  is 


=  i?sec^Z)»^C0.  4364  (H)  +  0.  5007Mi/i(H')/H*^  -  0.  0833/./  iH)/m 

Similarly,  by  adding  Eqs .  (16)  and  (21),  the  variance  of 
residual  wavefront  distortion  after  the  piston  and  tilt  are 
corrected  is 

Eirt—  ESm 

=  l-Jsec^Z)*^C0.  057/./  (/7)  +  0*  5007/.//,(^r)//7‘/*  -  0.  4513/./  im/m 


(22) 


(23) 


By  introducing  Eqs.  (22)  and  (23),  respectively,  in  Eq.  (1),  the 
expressions  for  dgpjj  after  the  piston  is  corrected,  and  dQpifp  after 
the  piston  and  tilt  are  corrected  can  be  derived  as  follows: 


^  i5-»/Scos‘^*^C0.  436/./  iH)  +  0.  5007/.//,(7f)/if’'*  -  0.  0833/./ 
iffft  —  ij"*/*cos’^*j>C0.  057/./  (H')  +  0.  5007 —  0.  4513/./ 


(24) 


(25) 


The  expression  for  dg  that  includes  the  piston  and  tilt  can  be 
derived  through  making  a  slight  change  in  Eq.  (6),  i.e.,  by 
expressing  part  of  atmospheric  turbulence  moments  as  omnidistance 
atmospheric  turbulence  moments,  and  then  inserting  them  in 
Eq.  (1)  as  follows: 


4o  =  ir*^cos' 


;3/S^( 


0.5007/.5/,\-*^ 


(26) 


3.  Atmospheric  Turbulence  Models 


In  November  1992,  Zeng  Zhongyong  et  al.  from  the  Anhui 
Institute  of  Fine  Mechanics  recorded  measurements  of  the  height 
distribution  of  Cp(h)  by  using  a  balloon-borne  exploring  meter 
in  Kunming  and  Beijing  and  obtained  huge  amounts  of  data.  We 
analyzed  the  measurement  data  from  Kunming  and  fitted  them  with 
the  Hufnagel  model.  The  Hufnagel  model  is  as  follows: 
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(27) 


Cj(i)  =«.i-exp(  -  i-)  +  «^xp(  -  Aj 

+  fl^xp(  —  ^  j  (m-*^) 

where  aj,  02,  Oj,  Hj,  H2,  Hj  and  a  are  all  fitted  constants;  h  is 
the  height  above  the  ground,  and  is  given  in  units  of  m.  Based 
on  the  time  of  measurement,  three  models  were  built,  namely 
morning,  night,  and  average.  The  fitted  parameters  are  shown  in 
Table  1,  of  which  the  data  from  Xinglong  were  derived  from 
reference  [7].  Fig.  2  shows  the  Kunming-average  measured  data, 
and  the  model-fitted  curves. 


,  2 

Fig.  2.  Cjj(h)  versus  height  at  Kunming  (in  Nov.  1992) 

The  squares:  the  measured  data;  solid  curve:  the  fitted  data 


4.  Calculations 


Taking  advantage  of  the  Kunming  atmospheric  turbulence  model 
and  HV5/7  model  introduced  in  the  previous  section,  we  calculated 
dQ.  The  calculations  were  based  on  both  Kunming  and  Xinglong 
average  turbulence  models;  in  addition,  for  comparison,  the  HV5/7 
[8]  model  was  also  calculated.  Fig.  3  (a),  (b) ,  and  (c)  show, 
respectively ,  the  calculations  associated  with  the  turbulence 


>8 
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Fig.  3.  Adaptive  optics  system  effective  diameter  dj  plotted 
as  a  function  of  laser  beacon  height  for  various  turbulence 
models  (a)  ,  (b)  and  (c) .  The  results  associated  with  the 

turbulence  model  Kunming-average,  Xinglong-average ,  HV5/7, 
respectively.  In  each  figure,  the  curves  correspond  to 
propagation  at  zenith  with  an  operating  wavelength  equal  to 
0.5  pm.  The  upper  curve  in  each  figure  is  for  the  case  in 
which  the  piston  and  the  tilt  are  corrected.  The  middle  is 
for  the  case  in  which  the  piston  is  corrected.  The  lower  is 
for  the  case  in  which  the  piston  and  tilt  are  included. 

model  Kunming-average,  Xinglong-average,  and  HV5/7.  During  the 
calculations,  the  zenith  ^  was  assumed  to  be  equal  to  zero  with 
an  operating  wavelength  A  equal  to  0.5|jm.  The  upper  curve  in 
each  figure,  derived  from  Eq.  (25),  is  for  the  case  in  which  the 
piston  and  tilt  are  corrected;  the  middle  curve,  derived  from 
Eq.  (24),  is  for  the  case  in  which  the  piston  is  corrected,  while 
the  lower  curve,  derived  from  Eq.  (26),  is  for  the  case  in  which 
the  piston  and  tilt  are  included  in  wavefront  distortion. 


It  can  be  seen  from  these  calculations  that  the  value  of  dj 
is  closely  related  to  the  distribution  of  Cu(h),  and  depends  on 
the  height  of  the  AGS  and  the  content  of  correction.  Let  us  take 
the  Kunming  turbulence  model  as  example;  if  the  height  of  the  AGS 
is  15km,  d|j  will  be  less  than  0.22m  before  correction,  and 
increase  to  approximately  0.58m  after  the  piston  and  tilt  are 
corrected.  This  means  that  with  an  0.6-m  aperture  telescope  and 
Rayleigh  scattering  at  the  height  15km  as  the  AGS,  the  variance 
of  residual  wavefront  distortion  can  decrease  from  5.3rad^  to 
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l.OSrad^  after  the  piston  and  tilt  are  corrected.  In  other 
words,  at  this  instant,  the  Strehl  ratio  of  the  system 
(SR=exp{-;E^}  can  attain  0.35  compared  with  the  0.0049  value  of 
the  Strehl  ratio  before  correction.  As  Close  and  McCarthy  Jr. 

[9]  confirmed,  when  SR~0.14,  a  star  image  that  approaches  the 
diffraction  limit  can  be  derived,  whose  half-peak  width  is  0.19'^. 
Hence,  it  is  possible  for  the  foregoing  system  to  derive  a  high- 
quality  star  image  in  Kunming. 

Finally,  it  is  to  be  noted  that  the  above  calculations  were 
conducted  with  reference  to  visible  light.  Since  dQ  is  directly 
proportional  to  the  calculations  are  also  applicable  for 

other  wavelength,  if  only  a  ratio  coefficient  is  multiplied. 
Obviously,  at  wavelength  farther  than  infrared  wavelength,  the 
effects  of  the  AGS  correction  will  be  even  more  remarkable. 
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